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Nature of Versatile Chemisorption on TiC(lll) and TiN(lll) Surfaces 
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Density-functional calculations on the polar TiX(lll) (X = C, N) surfaces show (i) for clean 
surfaces, strong Ti3d-derived surface resonances (SR's) at the Fermi level and X2p-derived SR's 
deep in the upper valence band and (ii) for adatoms in periods 1-3, pyramidic trends in atomic 
adsorption energies, peaking at oxygen (9 eV). A concerted-coupling model, where adatom states 
couple to both kinds of SR's in a concerted way, describes the adsorption. The chemisorption 
versatility and the general nature of the model indicate ramifications and predictive abilities in, 
e.g., growth and catalysis. 

PACS numbers: 68.43.-h, 73.43.Cd, 73.20.At, 71.20.Be 



In surface processes, like catalysis and film growth, 
and in materials design [through, e.g., chemical- (CVD) 
and physical- vapor deposition (PVD), or molecular-beam 
epitaxy (MBE)], initial adsorption is of key importance. 
Transition- metal carbides and nitrides (TMX; X = C, 
N) []| are widely used as substrates, coatings, and inter- 
layers in, e.g., cutting tools, electronics, catalysis, and 
biomaterials, and as components in the "MAX" phases 
y • A study of the nature of adsorption on TMX's is thus 
highly motivated. Generally, TiC and TiN are considered 
to be models for these materials. 
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This Letter reports on an extraordinarily rich ad- 
sorptive variation on TiX(lll) surfaces, found within 
a systematic density-functional-theory (DFT) study of 
CVD AlsOs/TiX coatings on cutting tools Q- Our 
trend calculations on adatoms in periods 1-3 of the 
table of elements show large variety in chemisorption 
strength and diffusion barriers (Fig. 1). To under- 
stand this, a detailed analysis of calculated bulk, sur- 
face, and adsorbate electronic structures is performed. 
The adsorption-energy (i?ads) trends are correlated with 
electron-structure trends: adsorbate-induced differences 
in surface densities of states [ADOS(i?)], real-space vi- 
sualizations of state-resolved DOS's [DOS(r,_E)], and 
valence-charge distributions, yielding a model for the 
nature of chemisorption on TiX(lll). The method re- 
lies on DFT calculations performed with the planewave- 
pseudopotential code dacapo Jj], with PW91 GGA 
exchange-correlation functional. Details are given in Ref. 
i 

Our analysis explains chemisorption on TiX(lll) as 
arising from the concerted coupling between adatom 
states and different types of TiX surface resonances 
(SR's) (Fig. 2): (i) a Ti-locahzed Fermi-level (Ep) SR 
(TiSR), which provides strong chemisorption through 
overlap with adatom frontier orbitals, giving well sepa- 
rated bonding (below the adatom level) and antibond- 
ing levels (above the TiSR and thus Ep) and (ii) X- 
localized SR's (XSR's) in the lower part of the surface 
upper valence band (UVB), which provide additional 
chemisorption, with strength depending on the energy 
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FIG. 1: Calculated atomic adsorption energies -Bads for fee 
site on TiX(lll) and for O in top-C site on TiC(OOl). Inset: 
energy differences between fee and bridge sites on TiC(lll). 



separation between XSR's and TiSR-modified adlevel. 
This concerted-coupling model (CCM) shows a deci- 
sive but earlier neglected role of subsurface X atoms for 
chemisorption. Its generality should make it relevant for 
other, similar, substrates. 

The NaCl-structured TiX's are covalent, ionic, and 
metallic |a|. The clean, polar, TiC(lll) surface, that 
is, prepared by flash heating under ultra-high vacuum 
(UHV) conditions, is Ti terminated (Fig. 2, inset) 0,1113 
and chemically very active, with dissociative adsorption 
of II2, O2, and methanol and a sticking for oxygen 100 
times stronger than on the stable, non-polar, and less 
reactive TiC(OOl) |l[l3,[l3- The TiC(lll) polarity is 
stabilized by a compensating surface charge provided by 
the partial filling of a strong TiSR ^i Ep II ES [I3- 
For TiN(lll), experiments are scarce, while theory indi- 
cates an environment-dependent termination [a, [XjI . To 
allow comparisons, Ti termination is here used for both 
(111) surfaces. 
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FIG. 2: (Color online). Schematic diagram of the concerted-coupling model (CCM), as illustrated by fee H on TiC(lll). 
Included are: bulk-TiC DOS(£'), clean-surface structure and D0S(J5), H/TiC ADOS(S) (c/. Fig. 3), and DOS(r,S) plots for 
(a-c) clean TiC(lll) and (d-f) H/TiC(lll), corresponding to the states indicated by the arrows [(b-c) are viewed perpendicular 
to the surface, (a) and (d-f) from above the surface, showing only the surface TiC bilayer]. 



For bulk TiX, our results confirm and extend earlier 
findings |3|. Strength, structure, and hardness result 
from directed iono-covalent Ti3fi-X2p bonds, with en- 
ergies in an X-dominated UVB (bulk TiC DOS in Fig. 
2). Corresponding antibonding Ti3d-X2p states lie in a 
Ti-dominated conduction band (CB), separated from the 
UVB by a pseudogap, a region of low but nonzero density 
of nonbonding Ti3d states. In TiC, Ep falls in the middle 
of the pseudogap (Fig. 2). In TiN, the extra electron per 
formula unit populates the CB, pushing UVB and CB 
to lower energies |5|, Ig and lowering the cohesion energy 
[14.8 (13.8) eV for TiC (TIN)]. Higher electronegativity of 
N makes TiN more ionic, yielding a higher X-localization 
of the TiN UVB H, and a higher Ti^X Bader [3 
charge transfer [1.5 (1.6) electrons for TiC (TIN)]. Our 
DOS(r,^)'s show that the bonding Ti3d-X2p states lie 
mainly in the middle-upper UVB region, while the lower 
UVB region is dominated by bonding X2p-X2p states 
[3, lial- Nonbonding X2s states form a separate lower 
valence band (LVB) at —9.5 eV. 

Fo^- TiX(lll), our DOS(^)'s (TiC in Fig. 2, TiN in 
Ref. 5) show significant Ti3d SR's around Ep, more filled 
on TiN due to its higher ionicity. DOS(r,^)'s [Fig. 2(a)] 



show that the TiSR's are strongly localized around the 
surface Ti atoms and extend toward the fee sites, i.e., 
what would be X sites in the bulk, there coupling with 
other surface Ti atoms. Thus, they are dangling bonds 
resulting from the cleavage of iono-covalent Ti-X states, 
which downshifts Ti-localized CB states and quenches X- 
localized states in the upper UVB region (Fig. 2). Cleav- 
age modifies also the lower half of the UVB. The sur- 
face DOS(-E)'s have three peaks below the main UVB 
peak, as opposed to one in the bulk (Fig. 2). DOS(r, £')'s 
show here a mixture of X-X bonding states, strongly lo- 
calized to the surface but also extending into the bulk 
[Figs. 2(b-c) illustrate this, with similar plots also for 
the other UVB peaks |a]]. Thus, the lower UVB region 
contains XSR's, causing peak-structure changes in the 
surface DOS(^)'s (XSR1-XSR4, see Fig. 2). 

Our adsorption calculations include full relaxation of 
all considered adatoms in fee, hep, and top sites on 
TiX(lll) and of oxygen in top-Ti and top-C sites on 
TiC(OOl). The calculated E^^s values jl^ show the fol- 
lowing trends (Fig. 1) [Sj: (i) much weaker O chemisorp- 
tion on TiC(OOl) than on TiC(lll); (n) for all adatoms 
on TiX(lll), preference for fee site, followed by hep; (iii) 



along both periods, pyramidic £^ads trends, peaking at 
group- VI elements (O and S); (iv) within each group, 
stronger bonds in period 2 than in period 3; (v) very close 
C and N values; (vi) very similar TiC and TiN values; 
(vii) a striking crossing of the TiC and TiN trends; and 
(viii) large variations in estimated fcc-hcp diffusion bar- 
riers (i.e., i^ads differences between fully-relaxed fcc/hcp 
sites and perpendicularly- relaxed bridge sites |18i|). 

A combined analysis of calculated ADOS(i<^)'s and 
DOS(r, i?)'s reveals the nature of chcmisorption. Key 
features are illustrated by the results for fee H on 
TiC(lll), showing (Fig. 2): (i) a strong negative Tid- 
derived double peak around Ep, corresponding to a TiSR 
depletion around the H adatom [Fig. 2(d)]; (ii) a strong 
positive Tid-derived peak at -1-1.1 eV; (iii) a strong neg- 
ative Cp-derived peak close to the substrate CSR3 peak 
(—3.9 eV); (iv) a broad Hs-derived positive region be- 
tween —6.1 and —4.1 eV (peaks at —5.5, —5.0, and —4.3 
eV), characterized by strong H~C bonding states [Fig. 
2(e)] and weaker H-Ti contributions; and (v) a positive 
region between —3.6 and —0.7 eV (peaks at —3.1, —2.3, 
and —1.4 eV), characterized only by H-Ti bonding states 
[Fig. 2(f)].^ 

This evidences a strong adatom-TiSR coupling that 
depletes the TiSR at Ep and yields antibonding states 
above Ep (at -1-1.1 eV) and bonding states below (pos- 
itive ADOS regions between —6.1 and —0.7 eV). Also, 
couphngs to CSR's are revealed, since (i) the strong 
broadening and splitting of the bonding adatom-TiSR 
states, (ii) the presence of negative peaks in the C- 
projected ADOS at the substrate CSR peak energies, and 
(iii) the presence of bonding adatom-C states in only the 
lower UVB region, successively stronger for lower-energy 
peaks, indicate that the positive ADOS peaks in-between 
the substrate UVB peaks correspond to combinations of 
bonding and antibonding adatom-CSR states (succes- 
sively more bonding for lower-energy peaks). Thus, the 
adatom couples in a concerted way to both TiSR and 
CSR's. This CCM is summed up in Fig. 2. 

Similar analyses for the other adatoms on TiC and 
TiN (Fig. 3 and Ref.0) corroborate the CCM. However, 
changes in adatom group number (III— > VII), period num- 
ber (2^3), and substrate (TiC^TiN) affect the calcu- 
lated electronic structures for adatoms in periods 2-3: 

(i) Increasing the adatom group number within each 
period downshifts the adatom energy level. The adlevel- 
TiSR energy matching is thus successively reduced. In- 
deed, the negative ADOS peaks at Ep decrease suc- 
cessively as group IV^VII (Fig. 3), indicating a weak- 
ening adlevel-TiSR coupling. On the other hand, the 
adlevel-XSR matching increases as group III^VI. Ac- 
cording to our ADOS's and DOS(r, -E)'s, this strengthens 
successively the adatom-XSR couplings, as exemplified 
by period-two adatoms on TiC(lll) [Fig. 3(a) and Ref. 
|5|: for B, the effective adatom p level lies at the upper 
UVB edge and is very sharp, indicating no coupling with 
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FIG. 3: (Color online). Calculated ADOS(£;) (difference in 
DOS for the surface TiX bilayer with and without adsorbate) 
for period-two adatoms in fee site on TiX(lll). Shaded re- 
gions mark the substrate UVB regions. Dashed vertical lines 
mark the energies of the substrate UVB peaks. 



UVB states; for C, the adatom p state lies in the upper 
UVB region, is slightly broadened, and forms a small, 
weakly bonding, adatom-C subpeak just below CSR4, 
indicating a weak adatom-UVB coupling; for N, the sit- 
uation is similar to H, the adlevel lying in the middle of 
the UVB and strongly broadening and splitting into sub- 
peaks of mixed bonding and antibonding adatom-C char- 
acter (more bonding at lower energies); for O, the adlevel 
lies in the CSR-dominated UVB region and forms strong 
adatom-C bonds. Similar results are obtained for period- 
2 adatoms on TiN(lll) [Fig. 3(b)] and period-3 adatoms 
onboth TiX(lll) surfaces 5]. Thus, the increasing group 
III^VI i?ads trend follows closely the increasing adlevel- 
XSR coupling trend, while very poor correlation is found 
with the adIevel-TiSR trend. The opposing adlevel-XSR 
and adIevel-TiSR trends can, however, explain the al- 
most constant E^ds value as C^N. 

(ii) Increasing the adatom period number within each 
group increases the adatom radius and upshifts the 
adlevel, due to core orthogonalization. The first ef- 
fect reduces adlevel-SR wavef unction overlaps. Our 
ADOS's show that this affects the adatom-TiSR cou- 
pling strength, decreasing the negative peaks at Ep J5|. 
The second effect changes the adlevel-SR energy match- 



ings. Our ADOS's and DOS(r, i?)'s show that this weak- 
ens the adatom-XSR couphng strengths, reducing the 
broadening and spHtting of the TiSR-modified adlevel 
13 • Thus, the decreasing period 2— >3 E'ads trend fol- 
lows closely these reductions in both adatom-TiSR and 
adatom-XSR coupling strengths. 

(Hi) Changing the substrate from TiC to TiN(lll) 
causes a downshift of all SR's and an increase of the 
clean-substrate TiSR filling. For a given adatom, the first 
effect reduces the adlevel-XSR matching, thus weakening 
the adatom-XSR couplings (Fig. 3). The second effect 
strengthens the adatom-TiSR coupling, as evidenced by 
the larger depletion of TiSR electrons on TiN (Fig. 3). 
For low adatom-group numbers, the adlevel-XSR match- 
ing is negligible on both surfaces and the adatom-TiSR 
coupling dominates, thus favoring adsorption on TiN. As 
the group number increases, the adatom-TiSR coupling 
weakens in favor of the adatom-XSR couplings, caus- 
ing successively stronger adsorption on TiC. Hence, the 
crossing of the TiC and TiN i?ads trends is well explained 
by the opposing trends in adlevel-TiSR and adlevel-XSR 
matchings as the group number increases. 

Thus, the i?ads trends are well correlated with electron- 
structure trends in which both types of SR's (TiSR and 
XSR's) are considered. This is further supported by: 

(iv) The weaker chemisorption at the end of each 
adatom period does not follow the energy-matching 
trends of (i) above. Yet, sharper ADOS p peaks for F 
and CI, compared to other adatoms with same adlevel- 
XSR energy matchings (Fig. 3 and Ref. y), indicate 
weaker adatom-XSR couplings for group- VII adatoms. 
The electronic structure shows a change in adatom- 
TiSR bond nature from iono-covalent to ionic as group 
VI^VII: smaller TiSR quenchings (Fig. 3), almost spher- 
ical adatom valence-electron distributions (Fig. 4), Bader 
charge transfers to the adatoms of 0.75-0.80 electrons, 
mostly from the nearest-neighbor Ti atoms, and mainly 
adatom-centered DOS(r, i5)'s at the adatom-peak ener- 
gies for F and CI. Thus, the weaker adatom-XSR cou- 
plings for F and CI are caused by their almost completely 
filled valence shells, which are essentially inert to interac- 
tion with the UVB. Also, for F/TiC(lll) the adlevel hes 
well below the substrate UVB, yielding further adatom- 
XSR weakening. These effects weaken the chemisorption 
as group VI— i-VII. 

(v) The much weaker chemisorption on TiC(OOl) evi- 
dences the importance of the adatom-TiSR coupling, as 
TiC (001) lacks TiSR at i;^ |li,|3. 

(vi) The stronger adsorption in fee than hep follows 
from the TiSR spatial extent [Fig. 2(a)], which favors 
adIevel-TiSR wavefunction overlap and thus adatom- 
TiSR coupling in fee. In turn, this causes a larger 
downshift of the TiSR-modified adlevel compared to hep, 
affecting the adlevel-XSR energy matching and thus 
also the adatom-XSR coupling strengths (as shown for 
0/TiC in Ref. 5). The largest such effect should be ob- 
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FIG. 4: (Color online). Isosurfaces of adsorption-induced dif- 
ferences in electronic charge densities (magenta, -f = increase; 
yellow, — = decrease; green = no change) for fee O and F on 
TiN(lll). The right panel shows a larger number of isosur- 
faces, illustrating also smaller absolute values of the charge- 
density difference, than the left panel. 



tained for adatoms with adlevel close to the XSR's. In- 
deed, our calculated fcc-hcp -Bads differences show pyra- 
midic trends that are very similar to those for the fee 
Bads values (with peaks, for TiC, at N and S |5|), show- 
ing again a very good correlation between calculated Bads 
values and adlevel-XSR matching trends. It can be ex- 
pected that diffusion-barrier variations (Fig. 1) are ex- 
plainable in a similar way. 

Thus, the CCM well describes the key features of 
the Bads trends on TiX. The large adsorptive variety 
is caused by the concerted action of adatom couplings 
to both surface TiSR and subsurface XSR's. The model 
should apply also to other systems with similar electronic 
structures, such as other TMX's, and thus have broad 
technological implications. 
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